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Abstract 

Using a sample of > 10^ J/^p — + fJ^'^ decays, Fermilab experiment 789 
has studied production of J/tp and V"' in 800 GeV proton-nucleon collisions. 
Differential cross sections and nuclear dependences have been measured for 
charmonium as well as for charm and beauty production. While charm and 
beauty production are consistent with perturbative QCD calculations, char- 
monium cross sections exceed the predictions of the color-singlet model by 
large factors, suggesting that additional mechanisms (such as color-octet pro- 
duction) may play important roles. Nuclear dependences of production cross 
sections may offer a new tool for the detailed understanding of charmonium 
production. 



1. Introduction 

Fermilab experiment 789 is a study of two-prong decays of beauty and charm 
which took data during the 1990/1 fixed-target run. E789 has published the most 
precise measurement of the nharm production cross section and its A dependence 
in 800 GeV proton coUisionstl as well as novel measurements at very forward xp 
the J/f/' production cross section and its A dependence using our beam dump (or a 
beryllium insert just upstream of the beam dump) as the target J3 More recent results 
include the only measurement to date of the cross section for beauty production 
in proton collisions at fixed-target energy,0 observed via the process tu^^ J/ip + 
X, J/ip ^'^d high-statistics studies of J ftp and ip' production.clla 
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2. Apparatus Description 

The E12S9 apparatus (shown schematically in Fig. ^ has been described in detail 
elsewhereQiJ; we summarize briefly here. It is based on the pre-existing E605/772 
spectrometer ,Q upgraded for this run by the addition of two silicon-microstrip vertex 
telescopes (Fig. ||), one above and one below the beam, replacement of the Station-1 
MWPCs with drift chambers, and a tenfold increase in data-recording capacity (to 
50MB/spill). Wire-like targets were used to localize the primary- interaction vertex 
in y and so that only the decay vertex need be reconstructed. 
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Figure 1: Plan (x-z) view of E789 apparatus. 
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Figure 2: Elevation view of E789 vertex telescopes; note that the target dimensions 
indicated correspond to the beauty running, and thinner targets were used for the 
charm running. 

The spectrometer features two large analysis magnets, SM12 and SMS, which 
deflect charged particles in opposite directions in y. A water-cooled copper beam 
dump suspended within SM12 absorbs non-interacting beam protons as well as sec- 
ondaries emitted within »i±20mr of the beam in y. Shielding within and around 
SM12 absorbs neutral secondaries. This geometry limits the pair acceptance to «1% 
but allows operation at high interaction rates. The vertex telescopes and the 23 
planes of scintillation-counter hodoscopes and drift chambers at Stations 1, 2, and 3 
measure the tracks of charged particles passing above or below the beam dump. The 
SM3 magnet serves to remeasure charged-particle momenta and thus to confirm the 
target origin of tracks, allowing the copious background of muons created within the 
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beam dump to be rejected. Particles are identified by electromagnetic and hadronic 
calorimeters, scintillation-hodoscope and proportional-tube muon detectors, and a 
ring-imaging Cherenkov counter. In most events only two oppositely-charged parti- 
cle tracks traversing the spectrometer are fully reconstructed, one passing through 
each vertex telescope. 

Data were taken separately in charm and beauty spectrometer settings. In the 
charm running we were able to operate the spectrometer at interaction rates up to 
5 MHz, using an on-line vertex trigger processorQ to reject >80% of hadron pairs 
from the target. The higher SM12 current used in the beauty setting allowed J /ip 
and data to be taken at a 50 MHz interaction rate, with no on-line vertex cut 
needed. 

3. Charm and Charmonium Production Cross Sections 

We have measured differential cross sections for charmEl and charmoniumO pro- 
duction. Charm data were taken using gold and beryllium targets at SM12 currents 
of 900 and 1000 amperes. Fig. ^ shows the observed hadron-pair mass distributions 
(under the K~'k~^ and 'k~ assumptions^) for the various E789 charm data sam- 
ples, using two different decay- vertex cuts. Clear D^{D^) signals stand out above 
the dihadron background. The looser vertex cuts (t/ctt > 7.2) were found to opti- 
mize the statistical significance of the signal and were used in the cross-section and 
A-dependence analyses. 




Figure 3: Left: dihadron mass spectra for the E789 charm data samples for two 
different decay-vertex cuts. Right: dimuon mass spectrum for E789 beauty data 
sample. 

§The RICH detector was not optimized for tlic D mass region and was not used in tliis analysis. 
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Our narrow range of acceptance in the longitudinal momentum of the pair (due 
to the beam-dump and shielding geometry) precludes a direct measurement of total 
cross section. We measured the D^/D^ differential cross section da/dxp = 58±3± 
7/ib/nucleon at (xp) — 0.03. This can he extrapolated over all xp (using the xp 
shape measured by previous. experimentsEl) to give a total / D'^ cross section a — 
17.7± 0.9 ± 3.4/ib/nucleonE Averaging with previous measurements using 800 GcV 
proton beamsB gives a{pN D° X) + a{pN -^D°X) = (20.aJi 3.5) ^b/nuclcon, 
consistent with next-to-leading-order (NLO) QCD predictionsEI within the broad 
range of theoretical uncertainty. 

Fig. ^ also shows the dimuon mass distribution from the beauty running, and 
Figs. ^ and ^ comcare our measured J ftp and tp' differential cross sections with 
QCD predictions.^ Here there is little uncertainty in the extrapolation over xp, 
with the xp shape well determined bothJpy these "open- aperture" data alone and in 
combinatiipn with our beam-dump data,a as indicated in Fig. ^. As at the Tevatron 
Collider ,EiJ J/tp and ^' production are substantially underestimated in the QCD 
calculation, with phenomenological "K factors" of 7 and 25 (respectively) required 
to give agreement in magnitude between data and theory. Note that the model 
calculation includes only contributions from color-singlet charmonium s tates and 
neglects possible contributions from color-octet charmonium componentsHiJIlj and 
from postulated states above DD threshold.Ej (These discrepancies have been the 
subject of much attention and are discussed further in Sec. 5.) We finda <7{p + N — *■ 
J/ip+X) = 442±2±88nb/nucleon_a|i(iYb+^ ^ -0'+^) = 75±5±22 nb/nucleon. 
Comparison with previous resultsta^ES (Fig. ^ shows that the J/ip total cross 
section and its excitation curve are by now well determined experimentally. 

4. Charm and Charmonium A Dependence 

Nuclear dependences of production cross sections can shed light on production 
mechanisms and thus are of intrinsic interest. In addition, suppression of char- 
monium production in nucleus- nucleus collisions has been proposed as a signature 
for quark-gluon-plasma formationJ23 so it is important to study processes respon- 
sible for charmonium suppression in proton-nucleus collisions, which might present 
a background to a quark-gluon-plasma signal. The production of heavy quarks is 
naively expected to depend linearly on the atomic weight (A) of the target nucleus, 
since the dominant X)CD mechanisms (gluon-gluon fusion ancLqg annihilation) in- 
volve hard partons.ESEI However, in Fermilab E772 we showedE3 that charmonium 
production in proton-nucleus collisions in fact has a complicated dependence on A, 
parametrized as da/dxp oc ^"(^j^)^ suggesting that other processes are at work in 
addition to those of perturbative QCD. 

The concentration of our charm sample in a narrow range of Feynman-x results in 
the most precise determination to date of the charm-production nuclear dependence 
at a point in xp, allowing a precise comparison to be made between open-charm and 
charmonium production. To augment the forward- a; A-dependence measurements 
of E772, we took additional J/ip data in E789, using a rotating wheel of beryllium, 
carbon, and tungsten targets placed 1.27 m downstream of the usual target position 
in order increase the acceptance near xp = 0. Fig. ^ shows the dimuon mass 
distributions thus obtained and the resulting A-dependence exponent a vs. xp for 
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Figure 4: Differential cross section per nuclcon vs. pt for J/?/; production com- 
pared with QCD prediction (note large phenomenological "if" factor by which the 
prediction has been multiplied). Dashed curves indicate indirect (via decay from 
higher-mass charmonium states) and direct contributions, and solid curve their sum. 

J/ip compared with that for D'^. We see that (at least at small xp) ct is significantly 
lower for charmonium than for charm: D*^ production depends linearly on AJxi = 
1.02 ± 0.03 ± 0.02 at (xp) = 0.03),EI while for the J/^j, a(0.03) = 0.89 ± 0.02.h13 

The increased nuclear suppression at low xp oi charmonium as compared to 
charm is consistent with models in which charmonium productioii is suppressed in 
nuclei due to dissociation by interaction with co-moving partons.Q Models in which 
the nuclear suppression of charmonium production is an initial-state effect {e.g. due 
to possible shadowing or nuclear modification of the gluon structure function) are 
disfavored, since they would predict similar nuclear dependences for .J/ip and D 
production. 

Effects on charmonium production at large xp due to intrinsic charmil (the. pres- 
ence of virtual cc pairs in the nucleon sea) and initial-state parton energy losgE3 have 
also been postulated. The predictions of the intrinsic-charm model have not been 
borne out by our data since they feature significantly larger (and more strongly 
A-dependent) cross-section contributions at the largest xp than are seenJj The 
qualitative trend we observe at large x p (nuclear suppression increasing with xj^l) is 
successfully accomodated in models which take account of parton energy losscJ in 
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Figure 5: Differential cross section per nuclcon vs. pt for ijj' production compared 
with QCD prediction (note large phenomenological "X" factor by which the pre- 
diction has been multiplied) . Also shown are phenomenological fits to the J/^ and 
ip' Pt shapes of the form (1 + (pt/po)^)~^ we find po = 3.00 ± 0.02 GeV/c for J/^p 
production and 3.60 ± 0.32 GeV/c for ■0'. 



traversing nuclear matter. Since gluons should interact more strongly with matter 
than quarks, the parton energy-loss model makes the intriguing prediction that at 
the highest Feynman-x (where qq annihilation dominates over gluon-gluon fusion), 
the nuclear suppression should decrease. Ea The apparent increase of a at ccp' ~ 0.8 
suggests that this may be occurring, but better statistics are needed for confirma- 
tion; these should be forthcoming from Fermilab E866, which is to take data during 
the 1996/7 fixed-target run. 

5. Discussion 

Production of charmEl and beautyi quarks is in reasonable agreement with per- 
turbative QCD calculations. On the other hand, production of charm and beauty 
quarkonia are observed at rates from one to two orders of magnitude higher than 
naively predicted. The experimental facts of enhanced quarkonium production in 
fixed-target experiments are not new, dating back to the mid-1970s when J/^ 
hadroproduction was first observed. At that time it was realized that lowest-order 
production of cc in a color-singlet state (via one intermediate virtual photon or three 
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Figure 6: Differential cross section per nucleus vs. xf for J /ip production in 800 GeV 
p-Au collisions along with phenomenological fit to the form (1 — jxi?!)"; we find 
n = 5.0 ± 0.2. (Note-that lowest-XF points in beam-dump sample are suspect due 
to large correctionsJij) 



virtual gluons) had diffiailty accounting for the large cross sections measured,c£l and 
the "color evaporation"L3 (or local-duality) mechanism achieved currency: the cc 
pair could be produced in a colored state and later emit a soft gluon to neutralize 
its color at little or no cost in probability. 

Recent advances in perturbative QCD have made the predictions of the color- 
singlet model computable with no free pa.ranieters,Ej allowing it to be definitively 
ruled out by both fixed-target and coUideiHO measurements. In the regime acces- 
sible to fixed-target experiments, questions may still remain as to the applicability 
of factorization and the role of intrinsic parton kt. But the collider data in the 
previously-inaccessible regime of pt 3> ?n (where perturbative calculations ought 
to be most trustworthy) have forced theorists to consider seriously additional non- 
perturbative mechanisms. As a result, the leading candidate models which hajsze 
emerged are the color-octet modeltilllj and an updated color-evaporation modelo 
These models, while less predictive than the color-singlet model, nevertheless make 
strong predictions (for example, that J/f/i's should be highly polarized at high pt), 
which can be tested in detail in upcoming experiments. In this connection we men- 
tion Fermilab E866, which should record > 10^ J /if) M^M^ decays (in closed 
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Figure 7: Energy dependencfi-of total J/?/' production cross section per nucleon; 
data from E789 and Refs. M-H. 



aperture) in the upcoming run, as well as the CO Charm projectJS in which a sam- 
ple of ^ 10^ J/tp decays could be accumulated, with (due to the open geometry) 
most final-state particles accompanying the J/^ also measured. Thus CO Charm 
holds the possibility of high-statistics measurements of Xc (as well as J/ip, and 
open-charm) production at fixed-target energy. 

Further data on A dependences could also be useful. Energy loss in nuclear 
matter may be a means to distinguish color-singlet and color-octet charmonium 
states, since the color-octet state has gluonic quantum numbers and may be strongly 
absorbed.L3 Distinguishing initial- and final-state A-dependence mechanisms calls 
for more data on the A dependence of open-charm production at large xp, another 
area in which CO Charm could contribute. 
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Figure 8: Left: dimuon mass spectra for E789 J /ip ^-dependence data samples. 
Right: The exponent a of the A dependence of the production cross sections for 
D'^(D^) and J/ip production vs. Feynman-x as measured in E789; also shown are 
J/ip resufts from E772. 
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